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Ghost inflation predicts almost scale-invariant primordial cosmological perturbations with rela- 
tively large non-Gaussianity. The bispectrum is known to have a large contribution at the wavenum- 
bers forming an equilateral triangle and the corresponding nonlinear parameter f'f^^^'' is typically of 
order O(IO^). In this paper we calculate trispectrum from ghost inflation and show that the corre- 
sponding nonlinear parameter t^l is typically of order 0(10*). We investigate the shape dependence 
of the trispectrum and see that it has some features different from DBI inflation. Therefore, our 
result may be useful as a template to distinguish ghost inflation from other models of inflation by 
future experiments. 

o . 
o 

■ I. INTRODUCTION 

^: 

' Almost scale-invariant primordial cosmological perturbations predicted by inflation fits observational data very 
well While this is certainly a great success of the general idea of inflation, there still remain many unanswered 
important questions about inflation. 

One of those important questions is how to distinguish different models of inflation observationally. There are many 
models of inflation which are consistent with observational data. It is often thought that tensor mode fluctuations [3| 
and non- Gaussianity 0, i, [Ml will be useful to distinguish some of them. Non-Gaussianity is the main subject of 
the present paper. 

' While single-fleld, simple slow-roll inflation predicts negligibly small non-Gaussianity [3], it is well known that there 
, are many ways to generate non-Gaussianities large enough to be detected by near future experiments. They can be 
categorized into three by epochs in which non-Gaussianities are generated: (i) super-horizon, (ii) horizon-crossing 
and (iii) sub-horizon epochs. Each of these three types has a bispectrum with characteristic dependence on shapes 
^ ■ of triangle formed by wave- vectors. The bispectrum for each type has a large contribution at (i) squeezed- (ii) 
\^ ] equilateral- [1-0] and (iii) folded-triangle Q, respectively. Among them, our interest in the present paper is on the 
I type (ii), in which large non-Gaussianity is typically due to higher derivative terms whose importance is enhanced by 
C ■ ' smallness of the sound speed. Concrete examples of this type includes k-inflation ^ , DBI inflation [1, [^ and ghost 
inflation [3, [H [H [H ■ 

While bispectnmr is the leading deviation from Gaussian statistics and thus a useful tool to distinguish some 
^ models of inflation from others, trispectrum can also provide additional information about inflation [9.-J4]. In this 

■ paper we shall calculate trispectrum from ghost inflation, hoping to find ways to distinguish ghost inflation from other 
inflationary models which predict similar bispectra. We shall find that the shape-dependence of trispectrum in ghost 
inflation shows some difference from that in DBI inflation [lol - [l^ . Therefore, the result of the present paper may be 
useful as a template to distinguish different models of inflation by future experiments. 

The rest of this paper is organized as follows. In Sec.[lT|we review the ghost inflation and show the powerspectrum 
d • and the bispectrum in this model. In Sec. IIIII we calculate the trispectrum in the ghost inflation. In Sec. IIVI is 
devoted to a summary of this paper and discussion. In Appendix |X] we give a brief review of the in- in formalism. In 
Appendix [B] we present the method to construct the interaction Hamiltonian. In Appendix [Cl we show the details of 
the calculations. 

II. REVIEW OF GHOST INFLATION 

In this section we briefly review the ghost inflation [7,] , which is an inflationary model in the ghost condensation [l5[ . 
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A. The Model 



The ghost condensation is the simplest Higgs phase for gravity and modifies gravity in the infrared. This can 
be realized if the derivative of a scalar field obtains a constant, timelike vacuum expectation value in a maximally 
symmetric spacetime, either Minkowski or de Sitter spacetime. By choosing the time coordinate properly, the vacuum 
expectation value of the scalar field can be written in the form 

{4>) = Mh, (1) 

where M is a constant. The background energy momentum tensor of the scalar field is identical to that of a cosmo- 
logical constant. 

Ghost inflation is a model of inflation driven by the scalar field </) responsible for ghost condensation. This situation 
can be realized, for instance, a la hybrid inflation. We introduce another fleld x, which is assumed to be massive for 
(j) smaller than some critical value (pc- At the critical point (j) — (j)c, x becomes massless and then becomes tachyonic 
for (f> > (j)c so that inflation ends as in hybrid inflation. An important point of this model is that the end of inflation 
depends only on the value of 0. 

In order to make predictions for the CMB anisotropy and compare with observational data, we need to consider 
perturbations around the background. Metric perturbations h^^, are deflned as 

ds^ = (-1 + htt)dt'^ + 2hudtdx' + [a^Sij + hij)dx'dx^, (2) 

where a = e^* . The scalar fleld is also decomposed into the background and perturbation as 

= Mh + TT. (3) 

The perturbation tt is the Nambu-Goldstone boson associated with ghost condensate, i.e. spontaneous breaking of the 
time reparametrization symmetry, and thus we call it ghostone [TBI] . Under the general inflnitesimal diffeomorphism, 
x'^ — >■ x'' + ^'^ , the perturbations tt and /i^i, are transformed as 

r2. 



■ij ^ hij + 2H^ta'^dij - di^j ~ djt 



h,j ^ h,j + 2H^ta'^d^j - d^^j ~ dj^^, (4) 



where the indices are raised and lowered by using the background metric. 

Let us now construct the quadratic action for tt in de Sitter background. We can do this in a systematic way fisj . 
We begin with the unitary gauge, i.e. tt = so that (j) = M'^t. This gauge condition still leaves residual gauge freedom 
corresponding to the time-dependent spatial diffeomorphism: 

t^t, x' ^ x' +C{t,x). (5) 

Under this transformation, the gravitational perturbations /i^^ are transformed as 

htt H> htt, hti H> hti - a^dt{a^'^ii), hij hij - di(,j - dj^i. (6) 

Since the quadratic action must be invariant under this transformation, it should be constructed from the following 
three terms: / dtdx^a^h^^, J dtdx^a^K^^ Kij and / dtdx^a^K"^, where 

K,, = i [a'dt{a-''h,,) - d.htj - djhu] , (7) 

K'l = a'^6^^5^^Kki and K = a-^S'^Kij. 

We can then obtain the quadratic action in general gauge by undoing the unitary gauge. This is achieved by 
performing the spontaneously broken diffeomorphism = M~^7r, under which 



htt ^ htt - 2M-^dtn, (8) 



+ Ar'^iHa'^S^jdtTT + d^djir). (9) 



After taking the decoupling limit M/Mpi — and thus dropping h^^i, in the quadratic action, we obtain the leading 
quadratic action for tt. 



5*2 = / dtdx^a^ 





, (10) 
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where ai and a2 are constants of order unity and we have normahzed tt. Here, = '^'■' Vj and = V'Vi. For 
/M^ <C 1, as we shall justify soon, we can drop the terms depending on HdtTr and thus the quadratic action is 
reduced to 



dtdx d 



3„3 



2AP 



(11) 



where a = ai + 02- Therefore, there is no (Vtt)^ term in the quadratic action for tt and the lowest spatial derivative 
term is a^''(V^7r)^. In particular, the dispersion relation for the ghostone tt is 

This implies that the fluctuation of tt oscillates in the regime fc^/M ^ Ha while it freezes for k'^/M <^ Ha. Moreover, 
because of the shift symmetry, interaction terms in the tt action always include derivatives and thus are suppressed 
outside the sound horizon, i.e. for k'^/M <^ Ha. Therefore, it is interactions during the epoch around the sound horizon 
crossing fc^/Af ^ Ha that essentially determines the bispectrum and higher-order correlation functions of cosmological 
perturbations. Near the sound horizon crossing, we have u! ~ (kjaf' jM ^ H and thus {k/aY ^ MH ^ H^ ^ Hw if 
H/M <^ 1. We shall see below that the COBE normalization for curvature perturbations fixes H/M to a small value 
([25|l . This is the reason why in (fTO|) . Hdtn can be neglected compared with a^^S/Wn. 

In order to see the leading cubic and quartic interactions, we should identify the scaling dimensions of time coordinate 
t, spatial coordinates and the ghostone tt. Suppose that energy E is scaled by a factor of s as £^ — > sE and that 
time t is scaled as i — )■ s^^t. Since we know from the quadratic effective action for tt that the dispersion relation is 
(X fc*, we have to scale fc as fc — >■ s^/'^k ov x as x ^ s~^l'^x. Then, by demanding that the action is invariant under 
scaling, the scaling dimension of tt is determined to be 1/4: tt — )■ s^^'^tt. With the use of these scaling dimensions, we 
can identify the leading cubic and quartic operators. 

When we introduced ([1]) as a consistent background, we implicitly assumed the shift symmetry, i.e. the invariance 
of the theory under a constant shift of the scalar field — > + const. Because of the shift symmetry, the fiuctuation 
TT appears always with derivatives. 

We further assume that the theory is invariant under the Z2 transformation cj) ^ —(j). This corresponds to the 
simultaneous change of the signs of t and tt: 



t 



-t 



and 



(13) 



For example, SittV^tt (without being multiplied by H) is forbidden by this symmetry. 

Let us now seek the leading cubic operator. The shift symmetry tells us that a cubic term has at least three 
derivatives. Because of the Z2-symmetry p^ . a cubic term should have one or three time derivatives. Combining 
these with the fact that the scaling dimension of spatial derivative is lower than that of time derivative, we conclude 
that the leading cubic operator is 



53 



dtdx^a^dtTT 



(Vtt) 



(14) 



where /3 is a constant of unity. This operator has the scaling dimension 1/4. Similarly, the leading quartic operator is 



8Af4 



dtdx a 



3 .(V7r)4 



(15) 



where 7 is a constant of order unity, and has the scaling dimension 1/2. 



B. Powerspectrum 

In this subsection we review the power spectrum of the curvature perturbation generated by ghost inflation. 
The scale M, which plays the role of ultraviolet cutoff of the effective field theory, is constrained by various 
observations. The strongest constraint is from the twinkling by lensing which bounds M < lOOGeV [l^ ^. Since 



It is known that the ghostone behaves Uke dark matter. If we suppose that the ghostone is responsible for all dark matter then the 
structure formation gives a lower bound on M: M > lOeV [TtII . 
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this scale is much lower than the Planck scale Mpi, decoupling limit is a good approximation. Therefore, in order to 
calculate the power spectrum of ghostone tt we can simply study the tt action without coupling to gravity. We shall 
later relate tt to the curvature perturbation C. 

The equation of motion for tt in the linearized level is 



2 ak^H'^if 



where Mk = aTTk, k is a comoving wavevector and k = |k|. Here, we have introduced the conformal time 



(16) 



ri = -- (17) 
a 

so that dt = adrj, and a prime represents derivative with respect to 77. 
We quantize u as usual 

Mk(i) = Wfe(t)ak + tUfe(t)aLk, (18) 

where ak and a'_j^ are annihilation and creation operators satisfying [ak,a]{j.,] = (27r)'^(5^(k — k'). Choosing Wk{t) so 
that it corresponds to the correct mode function in flat spacetime for very short wavelength, we have 

MV) = ^{~vy/'Hi]\iqr^'), (19) 

where 

(20) 

The observed fluctuations in the CMB are generated quantum mecanically and strechcd by the exponential expan- 
sion during the inflationary stage of the universe. Their wavelengths exceed the Hubble horizon and, at the end of 
the inflation, are much longer than the Hubble horizon scale. This means 

k\Ve\ « 1, (21) 

where k is the wavenumber of the mode of interest and rjf^ is the conformal time at the end of inflation. Therefore, it 
is a good approximation to take the limit rje — >■ 0. 

The power spectrum of the ghostone tt is obtained as 



" 27r2 



Wk 



" 7r(r(l/4)) W 



(22) 



The gauge invariant curvature perturbation C, 18| is related to tt as 



H H 



Therefore, we obtain the power spectrum of the primordial curvature perturbation as 



'^^ ~ ( M2 j " 7r(r(l/4))2a3/4 ) ■ (^4) 



The COBE normalization sets T'l''^ ~ 4.8 x 10"^, and thus 



^ ~ (1.6 X 10-3) X a^/io < 1. (25) 
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C. Bispectrum 



Let us now review the bispectrum of the curvature perturbation generated by ghost inflation. 

As in the previous subsection, we calculate a correlation function for tt in the decoupling limit and then relate 
the ghostone tt to the curvature perturbation C. For the calculation of the bispectrum, nonlinearity in the relation 
between tt and C, if any, should be taken into account. Fortunately, in our model of ghost inflation the simple linear 
relation (j23p is exact and there is no nonlinearity involved in the relation between tt and C, essentially because the 
background 4> is linear in t. Therefore, there are no additional nonlinear effects due to transformation of tt to C, and 
the bispectrum of ( is simply proportional to that of tt. 

As we have already seen in subsection III Al the leading cubic interaction in the tt action is 



/3 



2M2 



dtdx^a^dtTT^^ 



The corresponding operator in the Hamiltonian density is 



rLint,3 — O. r-^ OtTT ^ 

2M ^ 



(26) 



(27) 



Using the in-in formalism, we can calculate the tree-level bispectrum of ghostone tt, (tttttt), as 



dt' 



TTk^{t)TTk^{t)TTl^^{t), / dx^nint,3{t' 



(r(l/4))'a2 



j dr]r]-^F{i-OF (^^77^ F (^^77^ A:3(ki • ka) + symm. 



where 



(28) 



(29) 



Translating this into the bispectrum of C, we obtain 



(CkiCkaCka) = - 



_4\/2^3^ 

(r(i74)jV 



(2^3)<5M^k. n^." 



3? 



J d7jr]-^F{rj)F (^^^ F fc3(ki • ka) +symm 



(30) 



It is known that the bispectrum from ghost inflation has a large contribution at equilateral triangles formed by 

'Cqui 
NL 



three vectors ki, k2 and ks = — ki — k2 19]. The nonlinear parameter for equilateral configuration /^r^™' is defined 



as 



(Ck.Ck.Ck3)i.,..,=.3=fc - i^-mj:^^) ■ Ir^T' ■ 3 (^)" 

The value corresponding to the result ([30|) is (3] 



/X' =^ 85 • /3 • 



(31) 



(32) 



The WMAP 7-year data bounds f^^'^ [l 



as 



requil 
JNL 



26 ± 140(68% CL). 



(33) 



It is natural to take a ~ ^ ~ 1 as they are dimensionless parameters. In this case, the ghost inflation is not excluded 
by the WMAP 7-year constraint but can be detected in the near future by, say, the PLANCK satellite. 
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FIG. 1: The diagram for the contact term contribution. FIG. 2: The diagram for the scalar exchange contribution. 



III. TRISPECTRUM 



In this section we show the four-point function of late-time curvature perturbation ^ from ghost inflation. Details 
of calculations are presented in Appendix [C] 

As we have already seen in subsection III A[ the leading quartic interaction in the tt action is 



7 



8M4 



dtdx'^a' 



3 .(V^)' 



(34) 



As shown in Appendix [Bj the corresponding operator in the Hamiltonian density is 



T~LintA — 



8M4 



(35) 



where /? is the dimensionless coefficient of the leading cubic interaction. 

There are two contributions to the four-point function: a contact term contribution and a scalar exchange contri- 
bution. According to the in-in formalism (See Appendix E|. the leading contributions to the four-point function are 
written as 

(TTki (t)7rk2 (^TTka (i)7rk4 (^)> 



dt' 



d x'Tiint^it' , x'), TTfcj (t)7rfc2 (t)7rfe3 {t)Trki (t) 



dti 



dxl'Hint,3{ti,xi), / cixj'Hint, 3(^2, a;2), 7rki7rk2 7i'k3 7rk4 



(36) 



The first term in the right hand side of eg. ((36)) is called the contact term contribution and discribed in Fig. [l]as a 
four-point vertex. The second term is called the scalar exchange contribution whose diagram is shown in Fig. [2] This 
includes two three-point vertices. We calculate these contributions of the four-point function separately. 



A. contact term contribution 



Let us consider the contact term contribution to the four-point function. The contact term contribution to the 
four-point function is defined as 



(TTki (i)7i"k2 (i)7rk3 (^)7rk4 (i))cc = « / dt' 



1 



8M4 



^3^,(V7r(iV)r 
a{t') 



T^fcl (0^fe2 {t)T^k3 {t)T^ki (t) 



(37) 



which depends on 7. This parameter 7 does not affect the leading-order bispectrum and thus the four-point function 
provides independent information about the effective action. 

We define contact term contributions to trispectrums T^,cc('77 ki, k2, ka, k.^) and T^^cciv^ ki, k2, ka, k4) as 



{T^\.AtWAtWz{tWM) cc = (27r)353(ki -Hka -Hka -fk4)T^,cc(?7,ki,k2,k3,k4). 
(Cki(i)Ck.(i)Ck3(i)Ck4(i)>cc = (2^)3(53(ki -f k2 +k3 +k4)T^,cc(r/,ki,k2,k3,k4), 



(38) 
(39) 
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where r] is comoving time. As mentioned at the beginning of subsection III CI the simple hnear relation (|23p holds and 
there is no nonlinearity involved in the relation between tt and C. Therefore, we have the following simple relation 
between two trispectrums. 

r(;_cc(?7, ki,k2,k3,k4) = ^-^^ T^,cc(?7, ki,k2,k3,k4). (40) 

As shown in Appendix [Cl and presented as (|C9p . the contact term contribution to the trispectrum of curvaure 
perturbation is 



Tcccifl = 0, ki,k2,k3,k4) 



- / IT \ 12 / 
J I H \ I TT 



(91(729394) ((ki • k2)(k3 • k4) + (23 terms)) 



29 \M J Vr(i/4), 

{{-V^'Hi^iqW')) {{-Vr'Hi%{qW'))], (41) 

where 

9. ^ (42) 

^' 2M ^ ' 

The trispectrum depends on four 3-momenta k^ {i — 1,2,3,4) satisfying the constraint X)i ki = representing 
the momentum conservation. Assuming homogeneity and isotropy of the background, the independent variables are 
four amplitudes of the momenta, ki — |k.i|, and two angles between momemta, ki • k2/(fcifc2) and ki • k3/(A;iA;3). 
Although it is ideal to investigate the dependence of the trispectrum on all six parameters, it is somehow complicated. 
Fortunately, for the purpose of showing some differences between ghost inflation and other models of inflation such 
as DBI inflation, it is sufficient to investigate the dependence on a subset of six parameters. Following previous 
works [IQ . il2i] on the trispectrum from DBI inflation, we consider the equilateral case where 

ki — k2 — ~ — fc, (43) 

and the remaining independent variables are the two angles described above. As shown in (jC19p . the equilateral 
trispectrum of the curvaure perturbation is 

(1/2 \ 22/5 / \ 

zMf^) ifE^n^-^ (44) 

where 



, i=2,3A 



ki • k2 ks • k4 ki • k3 k2 • k4 ki • k4 k2 • k3 _ 

1,2 - 1,2 = '-^2, -L2 - 1,2 = '-'3' 1,2 " 1,2 = (.45) 

and Ci (i = 2, 3, 4) are constrained by the momentum conservation as 

1+ ^^ = 0- (46) 

1=2,3,4 

Fig E] shows the dependence of the contact term contribution to the equilateral trispectrum on the variables C2 and 
C3. In the most symmetric case where C2 = C3 = C4 = —1/3, the trispectrum has the smallest absolute value. 

Following Ref. [11], we define the nonlinecir parameter as 



nc . . ^ N rC,cc(kl,k2,k3,k4) -p|- 3 

r^, (ki , k2 , k3, k4) = (2,2^^(fc„))3 11 



NL 

4 



[(^1^2 + f^l^T) {^13 + ^14^) + (^1^4 + ^2^3) (^12*^ + ^13'^) + ('''1^3 + ^2^4) (^12"^ + ^14'^)] J 

(48) 
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where fcg is the wavelength coresponding to the present Hubble horizon size and kij = |ki +kj|. Since Vdk) is almost 
scale invariant, we do not distinguish between V({k) and Vc^{k[)) hereafter. In the equilateral case (|43)) . is 

We show the dependence of on the two angles C2 and C3 in Fig. |4l 
In the most symmetric case where C2 = C3 = C4 = —1/3, becomes 

r^L(fc,C2 = -l/3,C3 = -l/3,C4--l/3) = -1.007 xl04x A(^gM^!) ^ . (50) 



B. scalar exchange contribution 



In this subsection we present the scalar exchange contribution to the trispectrum. The scalar exchange contribution, 
corresponding to the diagram shown in Fig. [2J is defined as 



(7rki7rk2 7i'k3 7rk4; 



dU 



dti 



3l3a{ti) ( d 

dx^ — — 7rui,a;i 

^ 2M2 \dti ^ ' ' 



V7r(ii,xi) 



dr^- 

2AP \dt2' 



-^TT{t2,X2) \ (V7r(i2,a;2)) , Trk-.Trk^Trks'^ki 



(51) 



If the bispectrum ((30|) from ghost inflation is observed then it determines /3 and thus the scalar exchange contribution 
to the trispectrum. Then the full trispectrum from ghost inflation, which is the sum of this contribution and the contact 
term contribution described in the previous subsection, is fully specified by the other parameter 7 (or equivalently 7). 

Trispectrums TTr,se(?7, ki, k2, ks, k4) and T^^se(?7, ki, k2, ka, k4) are defined in the same manner as in the case of 
the contact term contribution. 



(vrki(t)7rk2 it)^^^,^{t)^Tl,^{t)) ^ 

(Ck,(i)Ck.(i)Ck3(i)Ck.(t)),e 
and are related to each other as 



= (2^)353(ki + k2 + kg + k4)T^,.e(ry, ki , ka , kg , k4 ) . 
(2^)353(ki + k2 + kg + k4)Tc,,e(?7, ki, k2, ks, k4). 



H 



7'C,se(?7, ki,k2,k3,k4) = ( — ) T^,se(?7, ki, k2, kg, k4). 



(52) 
(53) 



(54) 
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FIG. 5: This is the plot of the function A{C). 



FIG. 6: This is the plot of the function B{C). 



As in the calculations of other quantities, we take the limit i]e —>■ 0. Then, as calculated in Appendix [C] and presented 
in (|C23|) . we obtain 



7c,se(»7 = 0, ki,k2,k3,k4) 



M 



(27r)^(5^(ki+k2+k3 + k4; 



)) r-ri2 

— OO 'J —oo 



({ [(kx . k2)(fc? + kl + 2kx . k^) (gi^?)) fe^?)) {{-Vir^'H%Jq,,^l)) 

X [(k3 . U){kl + kl+ 2k3 . k4) {{~mr^'Hi%{q,r,l)) {{-mf^H^^q.^D) {{-m?^' H%^{q34vl)) 
~2{kl + ks ■ k,)kl {{-mf"H^Uq,^l)) {{-mf"H%^{q,r,l)) {{^^2^^ H^^^iqW^))] 

+23 permutational terms} + c.c.^ 
dr,, [(k, . k,){kj + kl + 2k, . k,){{-^,r/'Hi%{q,r,j)) {{~mf" H^^/\{q,^l)) ((~^i)^/^i?%ta2^?)) 



-2{kl + ki . k^)kl {(^Vlf/'Hi^iq^ril)) ((-^i)^/^i?%(<Z2^?)) ii-Vif^'Hi^iqW,))] 
dn2 [(ks • k4)(fc| + kj + 2ks ■ k^){{-V2f/'Hl%{q3vi)) {{-V2f/'Hl%{q4vl)) {.{-mf''H%{qW2)) 
-2{kl + ks • k,)kl {{-V2r/'Hi%{q,r,l)) {{-mr'H%^{q,^l)) {{-mf" H^^^^{qWM 



where qi is defined by eq. (j42|) and 



X {-r]w*^k^ (v)) {-Vwlk^ (77)) {-VWk:, iv)) {~VWki (v)) 

+23 permutational terms 

_ VZ^iJ(ki+kj)2 



2M 



, (55) 



(56) 



The scalar exchange contribution to the trispectrum Ti^^seiv = 0, ki, k2, ks, k4) depends on six parameters. As in 
the case of the contact term contribution, we consider the equilateral configurations by setting the amplitudes of all 
momenta to be the same (see eq. (j43p ). Then the remaining independent variables are ki •k2/ (fcifc2) and ki •k3/(fcifc3). 
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FIG. 7: The plot of the sum A{C2) + ^(Cs) + A{C4) 
as a function of C2 and C3 , where C4 = — 1 — C2 — C3 . 
The horizontal axis is rescaled by the value at C2 = 
C3 = -1/3. 



FIG. 8: The plot of the sum B{C2) + B{C-i) + B{Ca) 
as a function of C2 and C3 where C4 = — 1 — C2 — C3 . 
The horizontal axis is rescaled by the value at C2 = 
Cs = -1/3. 



As shown in eq. (jC30p . for equilateral configurations the scalar exchange contribution is reduced to 

1/2 \ 22/5 



r<;,,e(fc,C2,C3,C4) = 1.190 X 10- 



-16 



(n(fc)) 

4.8 X 10-5 



.8/5 



^ [A{C,) + B{C.,)]^ 

1=2,3,4 



(57) 



9/2 9/2 



nOO nOO 

A{C) = 72^(1 + Cf / dy2 / dyiyi^ 

Jo Jj^2 

X 1^(^3/4(2;?)) V_i/4 (2(1 + C)j/?) - if3/4(2/?)if-l/4(y?)if3/4 (2(1 + | 

X |c(x3/4(2/2))'^-i/4 (2(1 + C)yl) + ^^3/4(2/2)^-1/4(^2)^3/4 (2(1 + C)^^) | 
B{C) = {1 + Cf 

dy y^/^i C(i^3/4(y')) V_i/4 (2(1 + C)y2) _ K^/i{y^)K_^/^{y^)Ky^ (2(1 + 



(58) 



(59) 

where Ci {i = 2, 3, 4) are defined by eqs. (PS)) and constrained by momentum conservation as (HH). The scalar exchange 
contribution r(;^se(fc, (^2, C3, C4) is decomposed into six parts: those represented by A{Ci) and B{Ci) {i = 2,3,4). 
The functions A{C) and B{C) defined above are shown in Fig [5] and Fig [SI respectively. Fig [7] and Fig [5] show the 
sums A{C2) + A{C3) + A{C4) and B{C2) + B{C'i) + C{Ci), respectively, with C4 = -1 - C2 - C3 as functions of C2 
and C3. 

Fig. |9]is the plot of the total scalar exchange contribution Ti^^se{k, C2, C3, C4), which is the sum of the six parts. In 
the most symmetric case where C2 = C3 = Ci = — 1/3, the trispectrum has the largest value 



Tc,se{k, C2 — --, C3 — --, C4 



i). 2.561 X 10- 
3^ 



1/2- 



22/5 



4.8 X 10-5 



a- 



,8/5 



(60) 



This is always positive. Actually, Fig. [3] shows that the scalar exchange contribution to the equilateral trispectrum 
r^^se(fc, C2, C3, C4) is positive in the all parameter region. 



11 




FIG. 9: This plot shows the dependence of 
T^,se(fe, C2, C3, C4) on C2 and C3, where C4 = — 1 — 
C2 — C3. The horizontal axis is rescaled by the value 
at C2 = C3 = -1/3. 




FIG. 10: This plot shows the dependence of 
r^i(fc, C2, C3, C4) on C2 and C3, where C4 = -1 - 
C2 — C3. The horizontal axis is rescaled by the value 
at C2 = C3 = -1/3. 



FoUovifing Ref. [13|, we define the nonlinear parameter t^j^ as 

Jkr, k2, k3, k4) = ^2n^^^{k)r l\ ' 

X [(fc^A;2 + ^3^4) (^13^ + ^14'^) + {kfk^ + ^2^3) (^12'^ + ^13^) + (^1^3 + ^2^4) (^12"^ + ^14"^)] 

(61) 

As shown in eq. (jC32[) . t^j^ for the equilateral configurations is 



T^^i(fc,C2,C3,C4) ~ 8.945 X 10 



a«/^ U-8xl0-^; E.=2,3,4(l + C0-^/^ ■ ^ ' 

Fig. [TOl shows the plot of the shape of t^i^. 

In the most symmetric case where C2 = C3 = C4 = —1/3: 'Tnl becomes 

T^^Jfc, C2 - -1/3, C3 - -1/3, C, = -1/3) ^ 3.494 x 10^ x (^^^^^^) . 



(63) 



IV. SUMMARY AND DISCUSSION 



In this work we have calculated and investigated the trispectrum of curvature perturbation generated during ghost 
inflation. 

The analysis of scaling dimensions of operators makes it possible to identify the leading diagrams contributing 
to the trispectrum in ghost inflation. Actually, there are two leading-order contributions. One is represented by a 
diagram with one four-point vertex. This contribution is called contact term contribution. The other is represented 
by a diagram with two three-point vertices and called scalar exchange contribution. We have analyzed these two 
contributions separately. 

We have obtained general expressions for the two contributions as functions of six independent parameters. The 
six parameters are amplitudes of four 3-momenta and two angles between momenta. (Note that the sum of four 
3-momenta must vanish because of the momentum conservation.) 

In order to calculate the concrete values, we have focused on the equilateral case where all momenta has the same 
amplitude and where there remain two independent angular parameters as well as an overall amplitude of 3-momenta. 
Then we have calculated the non-linear parameters and t^j^ for the contact term contribution and the scalar 
exchange contribution, respectively. 
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In the case of local-type non-Gaussianity, it was forecasted that PLANCK will give the constraint \tnl\ ^ 560 [2l|. 
In the present paper we have shown that t^j^ and t^j^ are typically of order O(IO^). Therefore, the trispectrum from 
ghost inflation is probably detectable by PLANCK. Note, however, that the meaning of the non-linear parameter tjvl 
is different for different types of non-Gaussianities. (The same is true for the non-linear parameter f^^ of bispectra.) 
Therefore, as a future work, it is important to investigate detectability of the trispectrum predicted by ghost inflation 
in more detail. 

Now let us compare our results with the trispectrum from DBI inflation calculated in Refs. [H, E^l- The overall 
behaviors of trispectra are indeed similar. The trispectrum from ghost inflation has a peak at the equilateral config- 
urations, i.e. when all four 3-momenta have the same amplitude, because non-Gaussianity is mainly generated in the 
horizon-crossing epoch. This feature of trispectrum is shared with DBI inflation. Moreover, the dependence of the 
equilateral trispectrum on the angular variables C2 and C3 also has similarities. In both DBI inflation and ghost in- 
flation, the scalar exchange contribution has the maximum value at the most symmetric point C2 = C3 = C4 = — 1/3, 
and the absolute value of the contact term contribution becomes minimum at that point. The contact term contri- 
bution to the equilateral trispectrum has similar dependence on C2 and C3 in the two models of inflation. This can 
be easily seen by comparing Fig. [3] in the present paper and the right figure of Fig.l in Ref. p^. 

There are also some differences between the trispectrum from DBI inflation and that from ghost inflation. In DBI 
inflation, the value of the equilateral trispectrum is almost constant except for the edge region near the boundaries 
defined by = — 1 (z = 2, 3, 4). This feature can be seen in, e.g., the left figure of Fig.l of Ref. [l^. The trispectrum 
rapidly decreases near the boundaries Ci — —1 and the plateau looks like a triangle. On the other hand, as we can see 
in Figiniof the present paper, in ghost inflation the value of the equilateral trispectrum smoothly decreases towards 
the boundaries Ci = —1 {i = 2,3,4). As a result the shape of the plateau looks different from that in DBI inflation. 

If we look into actual values, we find an important difference. (Note that those figures mentioned above are 
normalized by the values at C2 = C3 = C4 = — 1/3 and thus do not tell the actual values of trispectra.) The sign 
and magnitude of the contact term contribution depend on the sign and magnitude of the quartic coupling constant. 
(This is in contrast to the scalar exchange contribution, whose sign does not depend on the sign of the cubic coupling 
constant.) In DBI inflation, the quartic coupling constant is determined by the sound speed. As a result, the equilateral 
trispectrum from DBI inflation has a positive value at the most symmetric configuration C2 = C3 = C4 = —1/3. 
On the other hand, in ghost inflation the dimensionless quartic coupling constant is an arbitrary parameter of order 
unity and thus the equilateral trispectrum at C2 = C3 = C4 = — 1/3 can be either positive or negative. Therefore, 
if equilateral-type non-Gaussianity is detected either by bispectrum or by trispectrum and if the negative equilateral 
trispectrum at the most symmetric point is observed, then it will support the ghost inflation scenario. 



Note added: While we were preparing the present paper, ref. |24| appeared on the arXiv. Before that, we had 
finished all calculations and presented our results at the workshop "The non-Gaussian universe" at Yukawa Institute 
for Theoretical Physics on March 25, 2010. The presentation file has been available from the workshop website at 
http://www2.yukawa.kyoto-u.ac.jp/~nlg/2010_3/prograni.htmsince March 27, 2010. 
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Appendix A: Review of in-in formalism 

In Heisenberg picture, the expectation value of an observable is 

{n\Q{t)\n), (Ai) 

where Q{t) is a Heisenberg operator corresponding to the observable and |r2) is the vacuum state in the initial time. 
A Heisenberg operator is sandwiched between initial vacuum states [22|]. 
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The time evolution of the Heisenberg operator is 

Q{t) = e'"'^*-'"^Q{to)e-'"<^'-'°\ (A2) 

where H is Hamiltonian. On the other hand, an operator in interaction picture evolves by the Hamiltonian Hfree of 
the free theory as 

Qj{t) = e*^/.ee(t-t0)Q(^^)g-iff/.„(t-t0)^ (^3) 

The relation between these operators is 

Q{t) = [Uit,to)]^ Qi{t)Uit,to), (A4) 

where 



= T <^ exp 



j' dt"Hi{t") 



(t > t'), (A5) 



Hi=H^Hfree, (A6) 

and T indicates time-ordered products. 

Let us now consider the time evolution of the vacuum state. First, we evolve the vacuum state 10)^^,^^ of -ff/ree by 
the total Hamiltonian H: 

n 

^ g-^Bo(to-*.„) 1^7) ^ _oo(l - ze)), (A7) 

where \n) (n = 1, 2, • • • ) represent excited cigenstatcs of H and 

H \n) = Eo \n) , H \n) = En \n) . (A8) 

Thus, 

|f])=C[/(to,-oo)|0)^_, (A9) 
where C is a constant factor. Similarly, we have 

=C*/,ee(0|[C/(^o,-Oo)]^ (AlO) 

By taking the inner product (r2|ri), we can determine CC* as 

1 = {n\n) = CC*. (All) 

Eventually, we have 

{n\Q{t)\n) = ^,ee(o|[[/(io,-oo)]tg(t)[/(to,-(x))|o)^,,, 

= free{0\ [U{t, -C^)]^ Q/(t)f/(t, -Oo) |0)^,,, 

= ^l^/ dtjv / dtM-l--- dhfree{0\[Hiih),[Hiit2),---[Hl{tN),Ql{t)]---]]\0)f^^^..iA12) 

J^—Q OO <y —OQ J —OO 

Appendix B: Interaction term of Hamiltonian 

If time derivatives of fields are included in interaction terms in a Lagrangian, the interaction terms in the Hamiltonian 
are not simply the minus the corresponding terms in the Lagrangian. The time evolution of a quantum state is driven 
by a Hamiltonian (See Appendix]^ but a theory is usually defined in terms of a Lagrangian. Thus, it is important to 
understand the difference between interaction terms in a Lagrangian and those in a Hamiltonian. In this appendix, we 
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consider this issue and obtain a perturbative expression for the difference. In Ref [ll|, this problem was investigated 
for a system with one particle. This appendix extends it to a multi particle system. 
We consider the case where the Lagrangian is of the following form: 



L(g^ (j^ t) = i free + Unt A",t), 



(Bl) 
(B2) 



where gab^ hab and Aab are independent of (7° (but can depend on time) and Lint{q'^,q'^,t) represents terms of third 
or higher order in q"- and q"". For simplicity, we consider the case where there is no constraint. In this case the inverse 
matrix of gab{t) exists and we denote it by g°'^{t). The canonical momenta are 



^ ^ - gab{t)q' + Aab{t)q' + ^^{q''A^{q\pb,t),t), 



(B3) 



and the Hamiltonian is 



H{t) = Paq''{q',Pe, t) - L{q%r{q',Pc, t),t) 

= Par{q',Pc,t) ~ ]^gab{t)q°'{q\VcM\'i\VcA + - L„,t(g", i), 

= ^3°' W(Pa - AUM){Vb - AM{t)q'') + habit)q\'' - L,„t(g^ r{q\pb, t),t) 

-lg''\t)^iq^,q'iq',Pd,t),t)^{q^,q'{q',Pd,t),t), 



(B4) 



where q'^{q'' ,ph,t) is defined by solving eq. (|B3|) with respect to 

Now, we quantize the theory. What we have to do is to promote the variables q° and Pa to the corresponding 
operators Q° and Pa- Then we have the standard commutation relations, 



[Q^n] = ^v, [Q'\Q'] = [Pa,Pb] = o. 

A complete set of quantum states is given by the eigenstates of or Pa which satisfy, respectively, 

Q^\q)^q\q), Pa\p)^p\p), 

and the inner products of them are 

{q'\q) = Y[Sq>^qa, (p'b) = nV„Pa, (qIp) ^Yi-^expiiq^-pa). 



(B5) 
(B6) 

(B7) 



In order to use the in-in formalism reviewed in Appendix [Xj we should adopt the interaction representation. The 
time-evolution equation of the operators Qf and Pi^a is a linear differential equation, where the subscript "J" implies 
that these operators are in the interaction representation. In the interaction representation, an operator Oj depends 
on time as 



Oiit) 



T 



Hfree{t')dt') 



0(to) 



T <^ exp(-i / Hfree{t')dt') 



to 



where 



Hfreeit) = ^Q^'^Pa - Aae{t)Q'){Pb - AdWQ'') + ^/^afeWQ^Q'. 



(B8) 



(B9) 



Thus, the time derivative of the operator 0/(t) is 



Oi{t) = i 



T exp(z / Hfree{t')dt') 



Hfree{t)0{ta) 



T<^exp(-i / Hfree{t')dt,) 



-i r|exp(i j Hfreeit )dt ) 
-l[Ol{t),Hfree,l{t)] 



0{to)Hfree{t) 



T < exp(-i / Hfree{t )dt 



(BIO) 
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The commutation relations (IBSP lead to the same commutation relations for and Pi a as 

[g?(t), = [Qm,Q]{t)] = [Pi At), PiAt)] = 0, 

and their time derivatives are 

PiAt) = -hah{t)Q\{t), 

Qm ^ 9''\t){PiAt) - Aumm)- 



(Bll) 

(B12) 
(B13) 



Since these equations are the same as those for free particles, we can easily represent Qj and Pi^a by using the creation 
and annihilation operators as in the case of free particles. 
The time evolution of a state represented as 



ij^m))i^H,^tAt)m))i' 



(B14) 



where 



Hint,l(t) — 



T'|exp(i^ Hfree{t')dt') 

H„AQ'i{t),PiAt)A, 



H,nt{Q%Pa,t) 



r<j exp(-Z / Hfree{t')dt') 
to 



(B15) 



and 



if»„t(Q", PaA= H{Q% Pa,t)- Hfree{Q'\Pa, t) 

= -L„AQ",r{Q\Pb),t) 

Solving eq. (jB3|) with respect to iteratively, we obtain 

81 

= - g''''^T^(Q^ Q^ t) + (more than third order in and Fa), 



(B16) 



(B17) 



where we have used eq. (IB13[) . Substituting this to eq. (jB16|) . we can see the difference between interaction terms in 
the Hamiltonian and those in the Lagrangian: 



1 f)T ar 

+ (more than fifth order in and Pa). 



(B18) 



Appendix C: Detailed calculations 

In this appendix we will show the concrete calculations of the trispectrum from ghost infiation. 

1. contact term contribution 



The contact term contribution to the four-point function can be written as 

7 f,,,iVnAA))' 



(TTki (O^Tka (OT^ks (07rk4 {t))cc = i I dt' 
Transforming variables as 



8M4 



X 



a{t') 



: TTki (^TTka (^TTka (t)i:i,^ (t) 



dt' = a{ri')dr], ttA) 



<i) 

Av'Y 



(CI) 



(C2) 



i I drj' ^ 
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we obtain 

(vTki (t)7rk2 (t)7rk3 (t)7rk4 (t)) 

r dry' / d^k[d^k'^d^k'^(k'^ ■ k'2)(-k; • k;, - • k;, - fc^,^) 

([(-^Vi(r/))(-7?'^ki(r/))(-V«k;,(r/'))(-'/"-k;-k-k;,(ry')), 

(-?7Mki (??)) (-Jy-^ks (f?)) (-»7*^k3 iv))i-VUk4. iv))] ) , (C3) 
where ui^iji) is written in terms of mode functions Wk{ri), w*_f.{ri) and operators Ok and a^j^ as 

Uk(?7) = Wfe(77)ak + wli.(77)al^, (C4) 

[ak,a],,] = (27r)3^3(k-k'). (C5) 
Substituting eq. (jC4l) into eq. (jC3|) . we have 
(TTki (i)7rk2 (i)7rk3 (0^k4 (i))cc 

7-8 /•>) 



.Tip 



/" d?7'(27r)353(ki + k2 + kg + k4){(ki • k2)(k3 • k4) + (23 terms)} 

^ — OO 

x{{-i^'w-k^{v')){-vw-k2{v')){-vw^k3{-n')){-ii'w^kA'n')) 

The contact term contribution to the trispectrum T^' ^^.(ki, k2, ks, k4) is 

r;',,,(ki,k2,k3,k4) = |^23?{*((ki • k2)(k3 • k4) + (23 terms)) 

i-VW-kiiv))i-W-k2iv))i-VW*_kArj)){~r]w*_kArj)) 

dri'{-j/w-ki{v')){~v'w-k2{v')){~v'w-k3{v')){-v'w-ki{v'))]- (C7) 



We are interested in the case where the physical length scales of all arguments ki {i = 1, 2, 3, 4) of the trispectrum 
are larger than the Hubble scale of de Sitter spacetime at the end of inflation, which means 

hrj < 1. (C8) 

Therefore, we approximate the trispectrum at the end of inflation as r.,^^*'(ki, k2, ks, k4). Then, we have 

(ki , k2 , ks , k4 ) = r;'.7," (ki , k2 , k3 , k4 ) 



^^^^4-g?|z((ki . k2)(k3 . k4) + (23 terms)) ^ (91929394)"'/' 23(r(l/4))-4 

ilTlj^' {hv^'H^Mv')) ((-.')'/^<i(.2.'^)) 

WM^ ( r(l"/4) ) (^i'?29394)"'/'' ((ki • k2)(k3 • k4) + (23 terms)) 

((-V)-V2i/«(93^''))((-r;')'/^ii^;i('?4r;'^))|, (C9) 
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where qi is defined in ea. p2)) and we have used 

V^-k, iv) - ~^{~vr^'Hi%im'), (CIO) 
^u;.,J,^„ = -z— (-j (Cll) 

Transforming variable as 

~V^V^x (<7 = <7i = 92 = 93 94), (C12) 
the contact term contribution to the equilateral trispectrum is 

r8 / ^ \4 / / \ \ ( M / / ^ \ 3/2 



where Ci is defined in cq.(|45|) and CiS are constrained by the momentum conservation as shown in ()46p . 

By using the relation (|23l) between <^ and tt, we obtain the contact term contribution to the trispectrum of curvature 
perturbation as 

(I, n n n \ ( \ ^-^^ecCfc, (^2, C3, C4) 

- (rf^)' (e/?) -'^{.f -(^--i;i<^^))'} . 

(C14) 

In order to tame the asymptotic behavior of the integrand at a; — > 00, we alter the integral route. The asymptotic 
expansion of Hunkel function is 



V ■nz 



2v + l 

I \ Z TT 



(z^oo). (C15) 



Thus, displacing the integral route from the real axis to an angle of 7r/4 does not change the value of integral. This 
change of the integral route is consistent with the ie prescription in the in-in formalism and picks up the correct 
vacuum. (See (jA7p .) We define this axis as y-axis, which means 

X = exp V- (CIB) 

Under this transformation of variable, we have 

X^/^ijW (X^) = ly3/2^^ (^2)_ (C17) 
^' ITT 

The expression (I24p for the power spectrum V(^{k) can be inverted to express H/M in terms of the power spectrum 
as 

^ ^ {j:{T{l/4)ra'/*Vcik)y^\ (C18) 

Hence, we obtain 

TcMk,C2,Cs,C,) = (^16y2^"/5(r(l/4))'/'5R|*^°°dyexp(Jz) {y^/' Ky,{y')] 

,11/5 7 



,1=2,3,4 



1/2 \ 22/5 



^2.215 X 10-" X ( (MIL 1 4_ I y Cf] k-^. (C19) 
\ 4.8 X 10-5 / I /-^ M ^ ' 



1=2,3,4 
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In the equilateral case where eq. (|43p is satisfied, becomes 

-1 



'^NL{k, C2, C3, C4) 



1=2,3,4 
1=2,3,4 



-3/2 



-3/2 



11/2 



E U- 



7 /-H^ 



87r2(r(l/4))4ai3/4 \m 



11/2 



r(i/4) 

(n(fc))-^5R|z (x^/^ijW (x^)) 



X]i=2,3,4 



E.=2,3,4 (1 + CO 



-3/2 ■ 



(C20) 



Substituting eq. (jCTSt into eq. ((C20)) . we obtain 

1 7 



'''NlQ^i ^2,6*3, C4) 



87r2(r(l/4))4 ai3/4 



(r(l/4))2a3/47'afc)) (T'clfc))-^?? I z r" dx (x3/^i?^;iC 



X]i=2,3,4 Cf 



I]i=2.3,4 (1 + Ci 



-3/2 ■ 



^(nW)^'/' Q^V^(r(l/4))^/^5i|z^°°d.(.V2i/«(x^))' 



Si=2,3,4 



.8/5 



E.=2,3,4(l + C.)"'/'' 



Si=2,3,4 Q 



E.=2,3,4 (1 + C^) 



-3/2 



-8/5 



-1.665 X 10^ X 



7 ( {Vdk)f'^ \ _ 
1^4.8 X 10-5 ) J2 



X]i=2,3,4 



i=2,3,4 



(1 + a) 



-3/2 ■ 



(C21) 
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2. scalar exchange contribution 



The scalar exchange contribution to the four point function is defined as 



dt2 / dti 



t I't2 

dt2 / dti < 



'2 2M2 \dU 



dxl [ -^T^{t2,X2) ) (V7r(t2,a;2))^ jTrkiTTkaTTkaTTk^ 



^^3_/3a(ti) f d 



2M2 \dti 



7r(ti,xi) (V7r(ii,xi))' 



'2 2M2 Wis' 



dxo ^,"'''^!^ f -r-7r(t2,a;2) ) (V7r(t2, 2:2))^ ) 7rki7rk2 7rk3 7rk4 



+ < TrkiTTkaTTkaTTki 



+ / dt2 ( dti 



dr 



3 Mil' 



1 2Af2 

TTki'I'ka 7rk3 7rk4 

Using eqs. (|C2I) . (|C4p and (|C5I) . we can transform this as 
(7rki7rk2 7rk3 7rk4)3g 



7r(ii,a;i) ) (V7r(ti, a;i))^ 



dx. 



3^a(t2 



X f dx^ ( ^ "(^2,a;2) \ / VM(7?2,a:2) y Mki(??) Mk2(r;) Mk3(^?) Mk4(??) ^ 

2 Vd772 0(772) / V a('72) / a(?7) 0(77) 0(77) 



+ / a{r]2)dri2 / 0(771)^771 ( / da;^ 



/ d u(7/i,xi)\ fVu{rii,xi) 



\drii a(77i) 



a(77i 



^^ Mki(?7) Mk2(??) Mk3(??) Mk4(?7) f ^^3 f d u{r]2,X2) \ f '^u{t]2, X2) 

a{ri) a{r)) 0(7/) 0(77) J \dr]2 0(772) / V aim) 



(C22) 
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-{4(ki.k2)(- 



+ {4(-fc2_ki 



+ 



-oo m 



dri2 



+{4(-fc2 - ki 



+ {4(ki-k2)(- 



+{4(-A;f - ki 



fl2 0-8 

: ^i^(27rWki + k2 + ks + k4) X 

— j — ^^{4(ki •k2)(k3 •k4)(-?7iu;_fe,(7?i))(-r?iw;_fc,(?7i)) (-r?iw;|ki+k2l('7i)) 

x(-r?2W-fe3(^2))(-r?2W-fe4(»?2)) (^^(-r?2t«*ki+k2|(??2))^ 

-'7wlfe, (??)) (-?7U'lfe2 (t?)) fc3 ('7)) (-'7w:lfc, (r?)) + 5terms} 

-{4(-fci - ki • k2)(k3 • ]s.4){-riiW-kA'ni)) (~'7i'*"-'=2 (^1))^ (-»7iW'|ki+k2| 

-'7U'-fci (??)) (-'7wlfe2 (11)) {-vwlks iv)) {-Vwlki (??)) + llterms} 
I - kg • k4) (-?7iw_fci (-77iw_fe2 (-»?i'«i'|ki+k2| (m))^ 

-?72W_fe3(r/2)) (^^(-r?2W-fe4(»?2))^ (-^2W*ki+k2|(»72)) 

-'7wlfe, (??)) {-riwlk2 iv)) {-VW-k3 iv)) {-Vwlk^ (??)) + llterms} 
k2)(-fc3-k3-k4)(-r?iw-/ci(r?i)) (-r^iW-fc^ (r/i))^ (-»?i^|ki+k2| (m)) 

-r?2«^_fe3(?72)) (^^(-??2W-fc4(%))) (-%w*ki+k2|(%)) 
r?w* (ry)) (-'ywlfc3 (77)) (-r?w* (r?)) + 23terms} + ex. 

{4(ki •k2)(k3 •k4)(-?7iw_fei(?7i))(-r/iM;_fc,(77i)) (^^(-mw'lki+kaKm)) 

-77wlfe^ (??)) (??)) (-?7Wfe3 (^7)) {-VWki iv)) 

-mwl^{r]2)){-V2'W*k^{m)) (^^(-%w^fki+k2|(%))) +5terms} 

k2)(k3 •k4)(-77iw-fci(??i))(-?7iw|ki+k2|(m)) (^(~^i^-'^2(»7i))^ 

-mwl^{V2)){-V2wl^{m)) (^^(-??2W^fki+k2|('72))) + llterms} 
I - k3 • k4) {-r]iw_k, iVi)) {-Viw-k2 iVi)) (-m«^|ki+k2| (?7i))^ 

-■nw'-k, iv)) {-vwlk2 iv)) {-v^ks (»?)) {-vwk^ iv)) 

-»72W^fe3(??2))(-ry2W|"ki+k2|(^2)) (^(~''2W^4(^2))^ + llterms} 

k2)(-fc3 -k3 • k4)(-r?i«;-fci(r?i))(-r?i«;|ki+k2|(m)) 

-Vwlki iv)) {-Vwlk2 iv)) {-VWk^ iv)) {-VWk4 (??)) 
-r?2W^3(?y2))(-?y2W^|ki+k2|('72)) (^(~^2W^4(^2))^ +23terms} 



^(27r)353(ki+k2 + k3 + k4)x 

-J ^({[(ki-k2)(-?7iu;-fei(r?i))(-77iw_fe,(r7i)) (^^(-r?iu>|ki+k2|('7i)) 

-2(fcf +ki •k2)(-r?iw-fei(77i))(-77iw|ki+k2|(?7i)) 

X [(kg • ]si4){-mW-hs{V2)){-V2W-kM) (-^U^fki+k^l (^2))^ 

-2(fc^ + k3-k4)(-772W-fc3('72)) (^^(-??2W_fcj772))^ (-??2«'|"ki+k2l('?2))] 

+23 permutational terms} + c.c.^ 



+ 



(/ ^{('*l''*^2)(-?7lW-fci(??l))(-m«^-fe2(m)) (^^(-??lW|ki+k2|(??l)) 



-2(A:^ + ki •k2)(-77iM;_fei(77i))(-?7iW|ki+k2|(»7i)) ( -^{-r]iW-kAvi)) )} 



/oo ^"^^ ("'^^^-'^3 (^2)) (-r/2W* fe, {V2)) (-r?2«^|ki+k2| (^2)) 



2(/c2 + k3 •k4)(-?72wlfe3(?72)) ( ;i^(-?72w!lfe,(7?2)) ) wfk.+k^ | (^2 )) } 



.'^??2 

X (ry)) (-r?w;* fe, (??)) (-??Wfc3 (ry)) {-r]Wk^ {v)) 



+23 permutational terms 



} 



(I)' (27r)353(ki + k2 + kg + k4)|- rfr?2 



'?2 

drji 



4M6 

({ [(ki . k2)(fc? + kl + 2ki . ^.,){{-n,f/'Hi%{q,r,!)) {{-Vif/'H^%{q2vl)) {{-Vif/'H%,{q,,r,!)) 
-2{kl + ki . k2)fc,2((-m)^/^/?Si(9i^?)) ((-m)^/^^f%(92^?)) {{-mf''H^^f,{q,,r,l))] 
X [(ks • k4)(fc3^ + kl + 2k3 • ^^4){{-mf''Hff,{q,4)) {{-V2f/'Hi%{q,nl)) {{-mf"H%,{qW2)) 
-2(fc| + k3 • U)kl{{-mf/'H^^/^{q,^l)) {{-mf''H%^{qW2)) {{-mf''Hf^^{qW.))] 
X (-?7wlfe^ (r?)) (77)) (-?7wlfe3 (r?)) (-ryw* (r/)) 

+23 permutational terms} + c.c. 



(1'^ rfryi [(kl . k2)(fc? + kl + 2ki . k2)((-^i)^/^i?(;i(giry?)) {{-Vi?"H^^f,{<l2vl)) {{-m?"H%,{q,,r^l)) 

-2{kl + \.^.\.^)kl{{-r,,f''H^i/,{q,^^^^^^ 

(^f^ dv2 [(k3 • k4)(fcl + kl+ 2k3 • \u){{-V2f/'Hi%{qsvl)) {{-V2f/'Hi%{q,nl)) {{-m?'' H%{qW2)) 

-2{kl + k3 ■ U)kl{{-mf"Hff^{q,r^l)) {{-mf" H%^{<l.vl)) ((-^2)^/^i/^Ji(<Z34%^))] ^ 

X (-??«;* fc, (7?)) (77)) (-7?u;fc3 (ry)) (-7?w;fc4 (??)) 



+23 permutational terms 



(C23) 
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where qi and qij are defined in eq. (|42p and eq. (|56p. respectively, and we have used the following identity: 

^ = 2/3ry {i~vf^'H%^iqr^')) {^ = 1, 2). (C24) 

Now let us take the limit — > and consider the equilateral case where the amplitudes of all momenta are the 
same. Then, we have 

TTT.se{k, C2, C3, C4) 

xxT[a{Hl%ixl)fH%i2il + C,)xl) - Hf^^{xl)H%,{xl)Hf^M^ + C,)xl)) + c. 



+ 

27V2 



' dx,xT{c,{H^^^,(xl)fH%^{2{l + - i?^)' (a:?)i/«/,(.T?)ijW (2(1 + 



(r(l/4))-^/?2a-i3/^i/Mii) {A{C.) + B{C.)l 



(C25) 



1=2,3,4 



where Ci, ^(C) and i?(C) are defined in eq. pS)) eq. (155)) and eq. ([55]) . respectively, and we have used eq. (IC17P and 

X-^/^H%^{x^) = ly-V2^_^/^(y2)^ (C26) 

x-^'^H%^{x^) = ^/2y-V2 (/_,/,(y2) _ ,^^/^(y2)) ^ (C27) 

(^2) ^ ^y3/2 (_/3/^(y2) ^ ,/_3/4(y2)^ _ (^28) 

Equations (jC28p and (IC27I) are special cases of the following relation: 

H'^\iz) = J,{iz) ~ iN,{iz) 

iCOSi/TT ^ ,. , i r ,■ \ 

= J^(iz) : Jv\^z) + J-u[iz) 

sm vir sm vir 



sm i^TT 
sin VTT 



\Ju{iz) + J-y{iz)) 



E 



(z/2) 



2n 



2/ ^n!r(i/ + n + 1) 

n— ^ ^ 



{z/2) 



2n 



21 ^n!r(-!/ + n + l) 

n— ^ ' , 



{-e^''"U{z) + I-,{z)) 



sm VTT 

The scalar exchange contribution to the four-point function of curvature perturbation is 

J-C,se[K, C-2, 03,04] = I ^ 



(C29) 



M4 



27^2 



(r(i/4))-^/3 



-4 /□2^-13/4 



' i=2,3,4 

215/2,6/5 (r(l/4))V5(^^(fc))lV5^^_9 ^ (A(Q) + S(a)) 



1=2,3,4 



1.190 X 10"^^ 



{ (n(fc)) 
I 4.8 X 10- 



1/2' 



22/5 



j,8/5 



k-^ ^ (A(Q) + 5(C0), 



(C30) 



1=2,3,4 
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where we have used eq. (jC18[ ) in the third equahty. The corresponding nonhnear parameter r^Lik, C2, C3, C4) is 



TNL{k, C2, C3, C4) 



V2 



E (1 + ^0" 

1=2,3,4 



■3/2 



T^^sejk, C2, C3, C4) 

(2^2pc(^o))3 



E (A(c.) + m)) 



i=2,3,4 



-3/2 



4=2,3.4 



7r7(r(l/4))4ai3/4 

Substituting cq. (IC18I) to this expression, we obtain 

24 /32 



94 R'^ / R \ ^^/^ 



i=2,3,4 



r^^^(fc,C2,C3,C4) ^ ^^j 4 J3/4 (^((r(l/4))^)a^/*7'c(fc))"^' (^W) 



E + 



-3/2 



1=2,3,4 



(C31) 



■(r(i/4)) 



i=2,3 
2/5 



a^/s V4.8 X 10 



-5 



E {A{a) + B{a)) 
=2,3,4 

E {A{a) + B{a)) 



E (1 + ^0^ 

i=2,3,4 



-3/2 



E i^ + cr' 



3/2 



1=2,3,4 



i=2,3,4 



E (1+^^) 



-3/2 



1=2,3,4 



(C32) 
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